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Abstract 

It is well known that attractive condensates do not posses a stable ground state in three di- 
mensions. The widely used Gross-Pitaevskii theory predicts the existence of metastable states 
up to some critical number of atoms. It is demonstrated here that fragmented metastable 

states exist for atom numbers well above . The fragments are strongly overlapping in space. 
The results are obtained and analyzed analytically as well as numerically. The implications are 
discussed. 
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Attractive Bose- Einstein condensates (BECs) have attracted much experimental 



and theoretical 
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The finding that the interaction between the atoms of a EEC can be varied systematically 
[2], allows one to prepare attractive condensates out of repulsive ones and to control them. 
A unique feature of attractive condensates in three dimensions is that they collapse, i.e., 
they do not posses a stable ground state as repulsive condensates do. At a given interaction 
strength, an attractive EEC can be formed, however, in a metastable state if the number of 
its atoms is smaller than some critical number N^r- 

The non-linear Schrodinger equation, also called Gross-Pitaevskii (GP) equation, is 



widely used to study condensates 
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221 1. This equation naturally predicts the collapse 

nrfp 

of attractive condensates and the appearance of metastable states, see, e.g., Refs. p, 18|, 19|, 
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22l |. The GP result can be derived using mean- field theory. Assuming a 
many-body wavefunction \E'gp(^i, ^^2, ■ ■ ■ , ^n) = 0o(^i)0o(^2) • • • <Po{^n), which possesses the 
structure of the wavefunction of N noninteracting indistinguishable bosons, one readily de- 
termines the GP equation by minimizing the energy expectation value Eqp = {"^gp \H\ "^gp) 
with respect to the orbital 0o(^)- Using, as usual, the Hamiltonian H = Hq + W, where 
Hq = h{i) describes bosons in an external trap and W the contact interaction between 



the bosons, the GP energy reads [2l|, |22| 



Egp = Nhoo + A( 



N{N -1) 



(1) 



where /iqo = (</'o|^|0o) and the interaction strength Aq is proportional to the s-wave scatter- 
ing length and takes on negative values for attractive condensates. 

The collapse and the existence of metastable states has been nicely discussed in the 

2l| as well as time- dependent Q, Q GP 



framework of the time-independent 
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equation which give identical results. Here, we follow the time-independent formulation. We 
consider without lost of generality a spherical harmonic trap and, as is common to all these 



approaches, we assume a normalized radial Gaussian function (foir) = (f 2 

[s^y/-^. For a = l, 



is the eigenfunction of 



for the orbital (f)o{r) = i^Y^'^'^o{r) with P 
a particle of mass m in the harmonic trap of frequency lj. Inserting this orbital in ([1]) one 
readily obtains 



Egp/ -hwN 
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+ (T^ + Ao 



A^- 1 



-a 



(2) 



where Aq = ^{'^Y^'^ f^'^^'^^;^- Since Aq is negative for an attractive EEC, this energy is 
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unbounded from below. For o" — > 0, the kinetic energy cr~^ grows, but cannot compensate 
the negative interaction term which scales as ct~^ , and the condensate collapses such that 
its density approaches a (5-function. Eq^{(t) exhibits a local minimum for < A*"^^ as is 
illustrated in Fig. [TJ To find the critical atom number , we set to zero the first and the 
second derivative of E^p^a) with respect to a, and obtain 
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in agreement with the result found using time- dependent approaches 



(3) 
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15| . The latter 



approaches also show that for a = r the r.h.s. of ([2]), i.e., Egp/N, plays the role of an 
effective potential for the motion of an atom at position r in the condensate. 

Within the GP theory, BECs cannot posses metastable states for > A"^^. The 



GP theory is a single-orbital mean-field theory in which all bosons reside in a sing. 



e or- 



bital 00 (r)- ^6 shall show below that within the multi-orbital mean- field theory [23|, |24| . 
which contains the GP theory as a special case, BECs do posses metastable states for 
A^ substantially larger than A'^^^. In this theory the many-body wavefunction is \1/ = 
^i) • • • 0o('"no)0i(^no+i) • • • 4'ii'^no+ni) ■ ■ ■ whcrc S IS the symmetrizing operator, and de- 
scribes the situation in which the A^ bosons are distributed among several orbitals: hq bosons 
reside in an orbital (poi^^, bosons in 0i(r*), and so on. Let us for transparency return to 
spherical traps and express the orbitals as (pimif) = ipi{r)Yim{0, ip). We consider the situation 
where the density p{r) is spherically symmetric, as is the case also for GP. Consequently, the 
occupation numbers riim corresponding to the (pim do not depend on m and may be denoted 
rii. The energy expectation value E = \H\ \E') now reads 



E = riohoQ + Sriih 
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\{po{r)f r'^dr + 







oo POO 

\ipi{r)\'^ r'^dr + GuqUi / \Lpo{r)f \Lpi{r)f r'^dr 
Jo 



(4) 
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where for simplicity we have considered the four orbitals spanned by / = 0, m = and 
/ = l,m = 0,±l. Of course, N = uq + 3ni. It is easily seen that this energy reduces to the 
GP energy Eqp in ([T]) if we put rii = 0. 

The many-body states described by the multi-orbital mean-field theory are fragmented 
states. In a fragmented state more than one eigenstate of the one-particle density operator is 



macroscopically occupied 



25|. In the GP theory the density reads PGp(r) — ^17^0(^)1^/(471") 



and as all atoms reside in a single orbital, this theory is unable to describe fragmented states. 
In the multi-orbital theory the density takes on the appearance 



(5) 



and the n/ describe the macroscopic occupations of the fragments. 
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29|. There, 



Fragmentation is by now well established in repulsive condensates 
as the repulsion between the atoms increases, the atoms minimize the energy by building 
fragments which avoid each other. At sufficientlystrong interaction, each atom may even 
reside in its own orbital in analogy to fermions [28 1. In the limit of infinite repulsion this 
behavior is known as the Tonks- Girardeau limit 30|. In princi ple, even the famous Mott- 



insulator state of atoms in wells of an optical lattice 3l|, l32|, |33l] can be viewed as a 



fragmented state, as each atom resides in its own orbital strongly localized in a different well. 
Indeed, the transition from superfluid to Mott-insulator states can be well understood within 



the multi-orbital mean-field theory 



3J]. In contrast to repulsive systems, fragmentation 



is rather counter intuitive in attractive condensates as the atoms attract each other and 



would like to be close to each other. There are, however, first hints in t 



le literature for 



fragmentation in attractive systems, the examples being atoms in a ring [29|, 135|, 136|, 137[ and 
in a symmetric one-dimensional double- well potential 23||. Interestingly, the fragmentation 
in these examples is dictated by symmetry and disappears when the symmetry of the trap 
is slightly distorted 23|, |35|]. The fragmentation in three dimensions discussed in the present 



work is not dictated by symmetry. We have performed the analogous computations for an 
anisotropic trap and came to the same conclusions as for the isotropic ones. 

We now return to (jlj) and compute the relevant quantities for the harmonic oscilla- 
tor. Each radial fragment orbital, in principle, contracts and collapses differently, and 
we give each ipi{r) a different scaling ai. The / = orbital is chosen as for the GP 
equation (with a = ctq) and for / = 1 we use the common respective radial harmonic 
'^i^'") ~ \f^'j7i('^^^^'^(^'^^~^^^^^'^^^^ ■ With this natural choice the energy expectation value 
now takes on the following explicit form 
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If one divides E by A^, one can immediately express the result for large N in terms of ctq, 
(Ti, AqA^ and the relative occupancy ^ only. 

At a given relative occupancy it is easy to compute the optimal values of ctq and o"i 
minimizing the energy E in (Q. An illustrative result for a cut along do of the energy E 
through its minimum is depicted in Fig. [T]for A^ > A"^^ and compared to the GP energy Eqp 
for the same Aq and A^. It is clearly seen that while GP theory predicts the collapse of the 
condensate, fragmentation does lead to the appearance of a metastable state. In addition, 
one notices that Eqp is larger than in ([6]) for a whole range of an ( see Fig. [1]). Interpreting 
E{cro = r) as an effective potential as done, e.g., in [sl, |lo|, H, 13, 15] for Eqp, this may have 



relevant implications also on the dynamics of an attractive condensate initially prepared as 
a repulsive condensate which naturally has a broad distribution (i.e., a large ao) and then 
switching the interaction to negative values. The thus created attractive condensate is likely 
to fragment as its momentarily lowest energy state is a fragmented one. 

The energy E is shown in Fig. [2] as a function of both ao and ai. This function exhibits 
a clear minimum. It is also seen that this minimum refers to a metastable state as E is 
unbounded from below for (Tq, cri 0. Numerically, it is straightforward to determine the 
critical number of atoms at which E{aQ, ai) in ([6]) ceases to support metastable states for a 
given fragmentation. In Fig. [3] this critical number N^r is drawn as a function of the relative 
occupancy Since N = nQ + Srii, the value of ^ can maximally be |. It is clearly seen in 
the figure that A"cr can be much large than A"^^. 

Although the dependence of E on ao and ai is rather simple, the analytical expression 
for A^cr is rather involved and not particularly informative. For this reason we computed 
the critical value N^^^"^ obtained from E m. ^ hy setting ctq = This can be done 
completely analogously to the calculation of A^^^ discussed above. One readily obtains for 
Ar> 1 

1 + 2^ 

iV-o=^i ^ ^GP ^JV (7) 

which gives || ~ 2.6 as an upper limit for the ratio N^^^'^^ / N^^ . The correct result for this 
ratio is very similar to the simple estimate ([7]) as can be seen in Fig. [31 

Having demonstrated that fragmentation leads to metastable states with A^ well above 
, we would like to give an ultimate numerical proof of this statement. To this end, 
we search fully numerically for a minimum of the energy functional (jl]) in the space of the 
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orbitals ipo{r) and v^i(r). These orbitals are thus determined self-consistently, similarly as 
done in the case of the single orbital in GP theory. Minimizing E[(po, ipi] leads to coupled 
equations for the orbitals which are solved numerically on a grid self-consistently. For 
details on the general multi-orbital equations and their numerical solution, we refer to the 



literature 
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34l |. In Fig. m we show the solutions (po{r) and (pi{r) for = 
l.TA*"^^, and compare them with those obtained using the "Gaussian model" ([6]) above. 
It is seen, in agreement with the model, that ipo and ipi simulate radial harmonics. The 
observed deviations are due to the full self-consistency of the numerical solution of the 
energy functional (jll). 

It is a rather very extensive task to explore the whole range of Ao,A^ and relative fragmen- 
tation using the full optimization of the energy functional E in i^. We have, therefore, 
performed calculations for several additional examples only. All the results obtained with 
this functional strongly support the findings of the "Gaussian model". In particular, we may 
conclude beyond doubt that fragmentation leads to metastable states of attractive conden- 
sates with particle numbers much larger than possible for unfragmented ones. In contrast 
to repulsive condensates where the fragmentation tends to separate the fragments in space, 
the fragmentation in attractive condensates is seen to be "on-top", i.e., the fragments prefer 
to be close to each other and to live in the same part of space. 

We would like to conclude with three brief remarks. We have numerical indications for 
a further increase of A^cr when allowing for an even higher order of fragmentation, i.e., by 
adding d-tjpe {I = 2) harmonics. The increase of A"cr is probably only moderate. We reiter- 
ate that the present findings are not restricted to spherical symmetry and harmonic traps. 
Fragmentation is expected to play a crucial role in the dynamics of attractive condensates, 
in particular of atom-rich condensates. 
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FIG. 1: (Color online) Energies per particle E/N as a function of the scaling parameter ctq. Shown 
are the GP energies for N < N^^ and N > N^^ (in black) and the energy of the fragmented 
condensate for N > (in red); see text for more details. is indicated on the respective 

curves. The relative fragmentation is ^ = 0.15. The scenario analyzed corresponds to *^Rb 
atoms in an isotropic harmonic trap of frequency uj = 2tt ■ 15Hz and s-wave scattering length 
of ag = — 7.14ao, where oq is Bohr radius, which are experimentally accessible {2]. For these 
parameters 

AtGP = 5000. The energy is expressed in units of huj and length in units of -o = -\/ — . 



Transforming to dimensionless units in which h = m = uj = 1, the corresponding interaction 
strength is Aq = 47ras/? = -0.001684. 
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FIG. 2: (Color online) The energy per particle E{aQ, ai)N as a function of Uo and ai (see Eq. ([6])) 
compared to the respective GP energy (see Eq. ([2])). = 6000. The quantities used are the same 
as in Fig. [TJ For these quantities the minimum of E is at ctq = 0.816 and ci = 0.811. 
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FIG. 3: (Color online) The critical number A'^cr obtained with E{aQ,ai) of Eq. ([6]) (black solid dots) 
and with the analytic result of Eq. ([7]) (red dashed curve) as a function of the relative fragmentation 
j^. The quantities used are the same as in Fig. [H N^^ is indicated by an arrow. 
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FIG. 4: (Color online) The orbitals (/5o(r) and ^i{r) obtained fully numerically by minimizing the 
energy functional ([H) for N = 1.7 (in black). The results are compared with those of the 
"Gaussian model" using Eq. ^ (in red). Here, ctq = 0.797 and ai = 0.788 at the minimum of the 
energy in Eq. The relative fragmentation is ^ = 0.25. The other quantities used are the same 
as in Fig. [TJ 
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